Helical spin textures in dipolar Bose-Einstein condensates 
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We numerically study elongated helical spin textures in ferromagnetic spin-1 Bose-Einstein con- 
densates subject to dipolar interparticle forces. Stationary states of the Gross-Pitaevskii equation 
are solved and analyzed for various values of the helical wave vector and dipolar coupling strength. 
We find two helical spin textures which differ by the nature of their topological defects. The spin 
structure hosting a pair of Mermin-Ho vortices with opposite mass flows and aligned spin currents 
is stabilized for a nonzero value of the helical wave vector. 
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I. INTRODUCTION 

Helical structures lie at the heart of several feats of in- 
novation. The Archimedean screw and spiral staircases 
are inventions dating back to ancient history. Simple 
bolts and spring coils stand as hallmarks of practicality 
in modern everyday life. Helical growth is a clear demon- 
stration that such structures also exist naturally without 
human intervention. On the other hand, helices are vital 
for life itself as we know it, the DNA polymer being a 
famous example of a double-helix structure [l| . 

The study of helical spin textures has recently drawn 
attention in the field of gaseous spinor Bose-Einstein con- 
densates (BECs). In a quantum degenerate gas of fer- 
romagnetic spin-1 87 Rb, a helical magnetization texture 
was observed to decay into small spin domains @- This 
effect was argued to result from the weak interatomic 
magnetic dipole forces in the system. Closely related to 
the experiment, the dynamical instability of an XY spi- 
ral state has been investigated theoretically ||. 

It has been predicted that a sufficiently strong mag- 
netic dipolar forces can spontaneously give rise to in- 
triguing spin textures in ferromagnetic condensates: the 
long-range dipolar potential stabilizes spin-vortex states 
in various geometries, as was demonstrated using a spin-1 
model [3, [1] and a classical spin approach [(| Even 
weak dipolar interactions can lead to dynamic forma- 
tion of a helical spin texture in a ferromagnetic spin-1 
BEC |8[. In the absence of external magnetic fields, the 
spin helix can appear as the ground-state texture in a 
suitable geometry and with strong enough dipolar inter- 
actions pj. 

In the present article, we investigate helical spin tex- 
tures in ferromagnetic dipolar condensates using a spin-1 
model. We find stationary helical states under the as- 
sumption that the system is infinitely long in the di- 
rection of the helix axis and find solutions hosting dif- 
ferent types of topological line defects, i.e., quantized 
vortices. Due to the symmetry of the order-parameter 



field, these line defects encircle the condensate in a heli- 
cal pattern. The resulting vortex structure resembles a 
pair of vortices excited in Kelvin modes in a stationary 
configuration. Direct experimental evidence of Kelvin 
waves has been observed in a spin-polarized BEC of 87 Rb 
by studying two transverse quandrupole modes of the 
atomic cloud Kelvin- wave excitations of quantized 
vortices have also been studied theoretically in a simi- 
lar setup H3-G3- Moreover, helical vortices in a two- 
component condensate have been investigated [liij . 

Although we focus on a spin-1 BEC, the study can 
also shed light on phenomena in strongly dipolar sys- 
tems with more complex order parameters. The energy- 
minimizing helical spin textures that we obtain are in 
qualitative agreement with those found earlier using the 
classical spin approach Q ■ Since the classical spin model 
is supposed to be accurate for ferromagnetic systems in 
the limit of large magnetic moments, it is plausible to ex- 
pect that similar structures exist independent of the value 
of the atomic spin number F. The most timely example 
of a condensate with larger F and significant dipolar in- 
teractions is the spin-3 gas of 52 Cr 1141 which has been 
produced by purely optical means [15j. The chromium 
atoms have magnetic moments of 6 [1b , whereas the maxi- 
mal atomic magnetic moment for an alkali-metal conden- 
sate is 1 hb ■ Recently, there has been progress in cooling 
and trapping vapors of Tm (4/xb) UMy Er (7/zb) [17j . 
and Dy (10 /^b) UM, the last having the largest atomic 
magnetic moment of all known elements. 



II. THEORY 



In this work, we study ferromagnetic spin-1 BECs, 
such as 87 Rb, using a zero-temperature mean-field model. 
In addition to the density-density and spin-spin inter- 
atomic forces, we also include long-range dipolar interac- 
tions in the model. 
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Stationary states of the system are solutions to the 
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time-independent Gross-Pitaevskii (GP) equation 
fctf (r) + ffn(r)*(r) + g s J2 M a (r)S a V(r) (1) 

a 

+ 5d ^ J dV^(r-r')M a (r')S^(r) = 0, 



a,/3 ' 



where ^ = (V>i, V'o, V>-i) T is a three-component spinor 
order parameter, h = — ft 2 V 2 /2m + Vt rap (r) — \x is the 
single-particle Hamiltonian, and $ a denotes the ath com- 
ponent of the dimensionless F — 1 spin operator whose 
spin-space expectation value gives the ath component of 
magnetization, M a (v) = * t (r)»S' a ^(r). The density of 
particles is given by n(r) = '0fe( r )V ; fc( r )i where the 
components (ipi, i/jq, ip-i) are the projections of 'f onto 
the eigenbasis of S z - The total number of particles per 
unit length, J d 2 rn(r) — N, is controlled through the 
chemical potential fi acting as a Lagrange multiplier. 

The coupling constants g, g s and gd measure the 
strengths of the local density-density, local spin-spin, 
and non-local magnetic dipole-dipole interactions, re- 
spectively. The first two are related to the scatter- 
ing lengths clq and ai into spin channels with total 
spin and 2h through g — 4irh 2 (ao + 2a2)/3m and 
g s = 4:irh 2 (a,2 — ao)/3m. Throughout the work, we 
use g s — —0.01 g, which is roughly the coupling con- 
stant for 87 Rb |19l - l21f and a value previously used, e.g., 
in 0, Hlj]. The dipolar coupling constant is given by 
9d = Mo^e^f/^ 71 " with no, hb, and gp being the per- 
meability of vacuum, the Bohr magneton, and the Lande 
factor, respectively. Rather than fixing gd to some par- 
ticular value, e.g., gd ~ 10 — 3 <7 as for 87 Rb, we present re- 
sults for various interaction strengths in order to empha- 
size the role of dipolar effects. Moreover, by using greater 
values of g&/g, the study should also provide useful infor- 
mation for systems subject to strong dipolar forces, such 
as gases of 52 Cr [2a |. In experiments, the ratio g^/g may 
be controlled with an optical Feshbach resonance [2J], 
which has been demonstrated for 87 Rb (25j . 

The long-range dipolar interactions are characterized 
by the functions D a p(R) = (S af3 R 2 - 3R a Rp) / R 5 , 
where {R a } denote the components of the argument R 
and R = \/^2 a R 2 ,- The elements of the traceless sym- 
metric tensor D a p have a <i-wave symmetric form and are 
thus simply expressed in cylindrical coordinates. The in- 
teraction integral in Eq. ([T]) can be viewed as an effective 
potential for the order parameter *S> at r arising from the 
magnetization throughout the system. 

In the present work, we concentrate on elongated he- 
lical solutions to Eq. (J]}, similar to the helical textures 
studied recently in cigar-shaped systems using a classi- 
cal spin approximation [7]. For simplicity, we assume an 
infinitely long system in the axial direction, confined ra- 
dially by a harmonic potential Vtmp = ^rnui 2 r 2 , where 
r = \/ x 2 + y 2 and u r is the radial trapping frequency. 
The axial symmetry of the trap allows a well-defined wave 
vector k for the helical texture. The aim is to fix k and 



calculate the energy-minimizing texture for a given ra- 
dial plane, say, for z = 0, subject to the condition that 
the planar texture is mapped along the axial direction 
according to the helical structure. Therefore, we write 
the order parameter in the general form 



V(r,(p,z) =e iq * z e iKzS *y(r,ip- 



KZ, 0) 



(2) 



expressed in cylindrical coordinates (r, tp, z). The second 
exponential factor describes a spin-rotation of angle kz 
about the z-axis in the clockwise direction. Due to the 
non-linear terms in Eq. ([1]), the order parameter feels 
an effective potential with a period of 2t:/k in the axial 
direction, and thus the Ansatz is written in the form of 
a Bloch-wave including the first exponential factor. 

By substituting Eq. ^ into Eq. (TTJ) and setting z = 0, 
we obtain a GP equation reduced to polar coordinates 

h'V(r, <p) + gn'{r, <p)*'(r, <p) + g s M a( r > V ) 

(3) 



where the primes denote that the quantities are evaluated 
at z — 0. The kinetic energy in the transformed single- 
particle operator hi is given by 



^2 r 



2m 



2m 



k L 



(4) 



where V 2 = d 2 + d r /r + d 2 /r 2 and L z = —id v . The 
functions Ip(r, if, k) arising from the dipolar interactions 
are evaluated in Appendix A. 

Only the kinetic-energy term in the energy functional 
depends explicitly on the axial wave vector q z , and hence 
minimization of the total energy with respect to q z yields 



K 

N 



d 2 r^ L z + S z ) * 



(5) 



Substitution back into Eq. Q reveals that for a fixed 
planar texture, the kinetic energy is a quadratic func- 
tion of the helix wave vector n. Moreover, the ki- 
netic energy becomes independent of k exactly when 
V&r(L z + Sz) = 0, implying that ^ must be an eigenstate 
of L z + S z . Denoting the (integer) eigenvalue by A, we 
find that in such a case the order parameter is of the form 
ipkiy, <Pt z) = e l ( A ~ fc ' ¥> -0/ c (r, 0, z), which is an integer-spin 
vortex state. This describes the general form of states 
for which the magnetization is invariant with respect to 
rotations about the z-axis, implying that the dipolar en- 
ergy and hence also the total energy are independent of 
n. 

Previous studies have shown that for very elongated 
condensates with finite dipolar interaction strengths, 
gd > 0, the spin-polarized textures with the magnetiza- 
tion lying along the long axis are the energetically favored 
ones [J, lZ| • Such states are of the general form discussed 
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FIG. 1: (Color online) A schematic illustration of a single pe- 
riod of a typical spin helix state. Here, p — 2n / k denotes the 
pitch of the helix whose wave vector is k. The solid curves de- 
pict the helical trajectories of two points on the z — plane, 
and the arrows point along the local direction of magnetiza- 
tion. If the state hosts vortices, their cores will also follow 
similar trajectories. 

above and hence do not depend on n. In this study, how- 
ever, we are interested in helical textures which are found 
by requiring the symmetry conditions [7[ 

M' x (r,<p + *) = -Mfr,<p). (6) 

Equations (|6]) are satisfied, e.g., when ipk(f, <p + 7r, z) = 
V'-fc( r > ft z )- The textures presented in the next section 
are solutions to Eq. ((1| subject to these symmetry condi- 
tions with respect to inversion about the z-axis. A typical 
spin helix is shown schematically in Fig. [TJ 

III. RESULTS 

We have solved the reduced GP equation, Eq. ©, for 
various values of the helical wave vector n and the dipolar 
coupling constant <?d- In the numerical simulations, we 
choose the harmonic oscillator length a r = ^/hjraui r for 
the unit of length and tkJ r for the unit of energy. For def- 
initeness, we fix the value of the dimensionless density- 
density coupling constant to g = gN/hu> r a^, = 10 3 , which 
corresponds roughly to 1.5 x 10 4 87 Rb atoms//zm in a trap 
with a r = 1 /im. 

One particular class of stationary solutions is illus- 
trated in Fig. [2] with <?d = 0.10 g. On the left-hand side, 
Figs.[2Ia)-[2][c), the state is shown in the axially homoge- 
neous case, k — 0, and on the right-hand side, Figs. EJ^d)— 
[5]Jf), with finite helical pitch, k = 0.25rad/a r . The left 
column in each subfigure shows the amplitudes of the 
order-parameter components \ipi\, \1po\1 an d IV'-iIj from 
up to down. The complex phases of the corresponding 
components are given in the right column in each sub- 
figure. For both states, the magnetization is pointing 
predominantly in the positive z direction in the upper 
half (y > 0) and in the negative z direction in the lower 
half (y < 0) of the xy plane. 
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FIG. 2: (Color online) Illustration of the MHV helix state for 
gd = 0.10g with n = in (a)-(c) and k = 0.25rad/a r in (d)- 
(f). For both states, the amplitudes of the order-parameter 
components (^i, ipo, ifi-i) are given in the left column and 
the complex phases in the right column. This class of states 
hosts a pair of MHVs with aligned spin currents and opposite 
mass flows. For each state separately, the color-map range in 
the left column is [0, max {n(r)} ] and in right column [0, 27r] . 
The field of view in each panel is [14 a r x 14 a r ] . 

The complex-phase plots in Figs. [2fa)-[2jc) reveal that 
the system is hosting two integer-spin vortices with phase 
windings (0,1,2) (upper half) and (—2,-1,0) (lower 
half) in the components (^i, ip®, ip-\ ), respectively. The 
ferromagnetic cores of the vortices are deformed into el- 
liptic shapes, which is also indicated by the separation 
of the phase singularities in the tp±i components, cf. 
Ref. [221 ]. The two vortices carry both spin- and mass 
currents: the spin currents flow in the same direction 
whereas the mass currents flow in the opposite directions, 
canceling the total mass current in the state. 

For finite helical wave vector, k > 0, the vortices 
move closer to the surface of the cloud, as illustrated 
in Figs. [Hd)-[2Jf). Not only do the cores of the vortex 
pair separate farther apart for larger k, but also the sep- 
aration between the singularities in the ip±i components 
increases. However, one should note that for larger k, the 
vortex lines are tilted steeper with respect to the axial 
direction. Due to the structure of the topological defects 
hosted by the order parameter, we refer to this state as 
the Mermin-Ho vortex (MHV) helix. This type of solu- 
tion is found to exist in the whole stability range of the 
system, < g d /g < 0.24. 

The MHV helix is the ground-state texture within the 
symmetry requirements of Eqs. ([6]) for all parameter val- 
ues considered in this work. However, there exists an 
interesting class of excited states which is illustrated in 
Fig. [3]for gd = 0.10 <? with n = and k — 1.0rad/a r on 
the left- and right-hand side, respectively. The complex- 
phase plots of Figs. Efa) and [UJc) reveal that the state 
hosts two spin vortices with the same phase windings 
(—1,0,1) in the components (-01, ipo, V'-i); respectively. 
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FIG. 3: (Color online) Illustration of the PCV helix state for 
gd ~ 0.10 g with k — in (a)-(c) and k = 1.0rad/a r in (d)- 
(f). For both states, the amplitudes of the order-parameter 
components (^i, tpo, ip-i) are given in the left column and 
the complex phases in the right column. This class of states 
hosts a pair of PCVs with aligned spin currents and vanishing 
mass flows. For each state separately, the color-map range in 
the left column is [0, max {n(r)} ] and in right column [0, 2ii\ . 
The field of view in each panel is [14 a r x 14 a r ] . 



These defects are spin vortices with polar (nonmagne- 
tized) core regions. Such vortices carry a spin current, 
whereas the mass current about the vortex line vanishes. 
The polar vortex cores within the ferromagnetic cloud are 
energetically analogous to air bubbles in water, hence in- 
creasing the total energy of the system through spin-spin 
interactions. The third singularity visible in Figs. [3ja) 
and [3fc) shows that the state hosts also a pair of frac- 
tional half-quantum vortices close to the surface of the 
cloud. However, these defects have very little effect on 
the texture because the related phase gradients lie in re- 
gions of nearly vanishing component amplitude \tp±i \ for 
all values of n. 

Solutions with similar vortex structures also exist for 
finite helical wave vectors, as depicted in Figs. [3]Jd)-[3]Jf) 
for k = 1.0rad/a r . Thus, we refer to this state as the 
polar-core vortex (PCV) helix. This class of stationary 
states is found within the dipolar stability range of the 
system, except for weak dipolar interaction strengths, 
9d/g 10~ 3 , likely due to numerical difficulties. 

If one considers the vortex pairs in Figs. [2] and [3] 
as single entities, the total phase windings in both 
states are (—2,0,2) in the order-parameter components 
{ipi,ipQ,ip-l), respectively. By tracing the order param- 
eter about the z-axis along a path close to the surface 
of the condensate, the expectation value of the spin ro- 
tates a total angle of 4ir about the axial direction. This 
suggests that both states represent doubly quantized spin 
vortices that have split into different kinds of singly quan- 
tized spin defects. The splitting of doubly quantized mass 
vortices in scalar condensates has been studied previously 
both experimentally and numerically 20 2!) . 



1 



(a) 


(b) 


1 — x — «- 




(<0 

/ mmm\ ^ * 

mm\\\\\\\\ i 1 1 


(d) 

, v v / ■ 

' < < • V V > 



FIG. 4: (Color online) Illustration of the spin textures for the 
states shown in Figs. [2] and [3] Panels (a) and (b) correspond 
to the MHV helix state with k = and k = 0.25rad/a r , 
respectively. Panels (c) and (d) refer to the PCV helix state 
with k = and k = 1.0rad/a r , respectively. The cones show 
the projection of magnetization onto the xy plane and the 
color indicates the local value of axial magnetization M' z (x, y). 
The field of view in each panel is [13 a r x 13 a r ]. 



The spin textures of the states depicted in Figs. [2] and 
[3] are shown in Fig. |4] Figures SJa) and B|b) illustrate 
the planar magnetization M'(r,tp) for the MHV helix 
and Figs. @[c) and [gd) for the PCV helix. The pro- 
jection of magnetization on the xy plane, "M.' (r, (p) — 
M' x (r,ip)e x + My(r,(f)e y , is depicted by the cones, with 
the length of each cone being proportional to the local 



magnitude ^ M^(r,ip) + My 2 (r,ip). The axial magneti- 
zation, M' z (r, ip), is shown by color. At the center of the 
system in the axially homogeneous textures (k — 0) in 
Figs. BJa) and |4jc), the magnetization is pointing per- 
pendicular to VM^ (r, (f) in the MHV and parallel to 
VM' z (r,ip) in the PCV state. For the MHV helix, this 
relative orientation is maintained for all values of the he- 
lical wave vector k. However, for the PCV helix, the 
relative orientation twists continuously with increasing k 
and finally locks into the symmetric configuration shown 
in Fig. |4|d) for k ~ 1 rad/a r . 

In Fig. [SJ the total energy per particle, E to t, is shown 
as a function of the helical wave vector k for the MHV 
and the PCV helices. As previously, the coupling con- 
stants have the values g = 1000 and <?d/.9 = 0.10 for 
both states. The MHV helix is the ground-state tex- 
ture, within the constraints of Eq. ©, for all values of k. 
Whereas E tot for PCV helix is minimized for the axially 
homogeneous state, k = 0, the MHV helix is stabilized 
for a finite wave vector K m ; n . This minimum persists 
in the total energy for all positive interaction strengths 
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FIG. 5: (Color online) (a) Total energy per particle, E tot , as a 
function of the helical wave vector k for the MHV helix (solid 
curve) and the PCV helix (dashed curve) with g = 1000 and 
gd/g — 0.10. The total energy of the MHV helix exhibits a 
local minimum with a finite wave vector K m ; n . (b) The energy- 
minimizing wave vector K m in for the MHV helix as a function 
of the dipolar interaction strength gd/g for two values of the 
density-density coupling constant, g = 1000 (solid curve) and 
g = 200 (dashed curve). 



g and <?d considered. However, the relative magnitude 
of the dip, [JSt o t(0) - .Etot(«nun)] /E tot (0), decreases for 
weaker dipolar interaction strengths. The curve E to t(K) 
is roughly a shifted parabola, closely resembling the re- 
sult that long-period helical structures in MnSi and FeGe 
can become stable due to ferromagnetic Dzyaloshinskii 
instability [30j . 

The energy-minimizing helical wave vector K m ; n is 
shown in Fig.[3Jb) as a function of the dipolar interaction 
strength gd/g for two values of the dimensionless coupling 
constant, g = 1000 (solid curve) and g — 200 (dashed 
curve). In the absence of dipolar interactions, the total 
energy is minimized by the axially homogeneous texture 
for which the kinetic energy is minimal. The value of K mm 
increases rapidly as a function of gd/g because the dipolar 
interactions, which overwhelm the kinetic energy already 
for gd/g ~ 0.01, favor a finite helical pitch. In general, 
Kmin increases for stronger dipolar couplings. However, 
the pitch of the energy-minimizing helix, p nlln — 2ir / 'K mm , 
is not determined directly by the dipole-dipole coherence 
length £dd = h/ y/2mgdn(0), where n(0) is the particle 
density at the trap center. In fact, the coherence length 
£dd shrinks for larger particle numbers, whereas p m i n typ- 
ically increases. This is likely due to the overall expansion 
of the condensate. 

As a measure of how large the axial magnetization 
M z (r) is on average, we define the integrated axial mag- 
netization per unit length as M z = J d 2 r\M' z (r, ip)\. Sim- 
ilarly, the integrated transversal magnetization is given 

by M r = J d 2 rJM^(r, ip) + M®(r, cp). These quanti- 
ties are plotted in Figs. [Ua)-[5Jd) as functions of the 
helical wave vector k for g — 1000 and four different 
values of the dipolar coupling constant gd/g. The solid 
curve corresponds to M r and the dashed curve to M z . 
As shown in Fig. [SJa), the axial magnetization vanishes 
when gd/g — 0. For k > 0.7rad/a r , the system becomes 
nonmagnetized because the additional kinetic energy due 




K [rad/a,.] 1.5 k [rad/a r ] 1.5 



FIG. 6: (Color online) Integrated transversal A4 r (solid curve) 
and axial M z (dashed curve) magnetization per unit length 
as a function of the helical wave vector k for the MHV helix. 
The strength of the dipolar coupling gd/g is indicated in each 
panel. The overall magnetization is reduced for larger values 
of k, and the axial magnetization tends to increase with gd/g- 



to finite k for magnetized states exceeds the energy gain 
from the spin-spin interaction, cf. Ref. 0. Already for 
weak dipolar coupling, gd/g = 10~ 3 , which is roughly 
the value for 87 Rb, the axial magnetization becomes sig- 
nificant, as shown in Fig. E^b). Also, for finite dipolar 
coupling, the value of k at which the system enters a 
polar state is increased. The axial component of magne- 
tization becomes larger for increasing gd, as illustrated 
in Figs.[6jc) and[6jd). 



IV. DISCUSSION 

In summary, we have studied helical spin textures in 
a spin-1 BEC subject to long-range dipolar interactions. 
The axially elongated system was assumed to be con- 
fined radially by a cylindrically symmetric harmonic po- 
tential. By using a helical Ansatz, we reduced the zero- 
temperature GP equation to a two-dimensional problem, 
in which the helical wave vector k appeared as a param- 
eter. This allowed us to investigate states with variable 
values of k. We found two classes of helical solutions 
which we refer to as the Mermin-Ho and the polar-core 
vortex helices, according to the structure of their topolog- 
ical defects. Whereas the total energy of the polar-core 
vortex helix is minimized for the axially homogeneous 
case, k = 0, the Mermin-Ho vortex helix is stabilized for 
a finite pitch. 
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The helical spin textures studied in this work are nat- 
urally most transparent in elongated systems, such as 
cigar-shaped BECs. One difficulty in observing them 
as stable configurations in a condensate subject to weak 
dipolar interactions is that the spins tend to align pre- 
dominantly parallel to the weak axis of the trap due to 
the head-to-tail attraction of the dipoles [3, El 0] ■ How- 
ever, we have performed additional three-dimensional 
simulations indicating that this problem can be overcome 
by placing the cigar-shaped system in a one-dimensional 
optical lattice potential. A strong enough lattice deforms 
the elongated condensate into a series of oblate clouds. 
Within each cloud, the preferred direction of magnetiza- 
tion lies perpendicular to the weak axis of the trap, and 
the relative orientation of magnetization between neigh- 
boring systems is determined by the long-range part of 
the interaction potential, i.e., by dipolar forces. 

In the three-dimensional simulations, we calculated the 
total energy of an axially spin-polarized state and a he- 
lical spin texture subject to the symmetry constraint in 
Eq. ^ as a function of the strength of the optical lat- 
tice potential. We chose the aspect ratio of the confin- 
ing harmonic potential such that uj z /uj r — 0.10, where 
lo z and uj r are the axial and radial trapping frequencies, 
respectively. The values of the dimensionless coupling 
constants were fixed to g — 5 x 10 4 , g s /g = — 10~ 2 , and 
<7d/<? = 10~ 3 , corresponding to 8 x 10 5 87 Rb atoms in a 
trap with u r = 2tt x 100 Hz. The optical lattice potential 
was of the form V op t — E s sin 2 (k s z), where the constant 
E s is the strength of the lattice potential, and the wave 
vector k s determining the distance between neighboring 
clouds was fixed to k s = n/10a r . The energy differ- 
ence between the helical and polarized states was found 
to decrease monotonously as a function of E s , reaching 
degeneracy at E s « 17 huj r . 
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Appendix A 

Here, we calculate the dipolar integrals 

Ip{r,tp,K) = Y, I d 3 r'D afj (r-r')\ z=0 M a (r') (Al) 

a J 

appearing in the reduced GP equation, Eq. ([3]). As shown 
below, the dipolar potential for a given magnetization 
can be efficiently evaluated by carrying out a series of 
one-dimensional Hankel transformations. 



The expectation values of the spin operators S a are 
readily evaluated for the Ansatz in Eq. © by applying 
the Hadamard lemma. Recalling that quantities evalu- 
ated in the z — plane are denoted by primes, we obtain 

M x (y) = cos(nz)M x (r, ip + kz) + sin(nz)M' v (r, tp + kz), 
M y (r) = — sm(Kz)M' x (r, ip + kz) + cos(Kz)M y (r, tp + kz), 
M z (T) = M' z (r,ip + Kz). (A2) 

The components of the planar magnetization are then 
expanded in polar Fourier series as 

Af;(r,^)=5>; l (r) e "T (A3) 

■n 

where 5 ™(r) = / dcpM^r, ip)e~ mv /2k = [ 5 ""(r)]*. For 
the helical textures, the components M' are in general 
slowly varying functions of the azimuthal angle <p, and 
hence they are accurately approximated by only a few 
terms in the expansion. The components of magnetiza- 
tion can be written in the form 

M a (r) = J2 C s ai g^{r)e^ KZ+n ^ + ^\ (A4) 

s,n,7 

where s G {1,0,-1} and the coefficients arising from 
the trigonometric functions in Eqs. (|A2[) are given by 
@xx = = s 2 /2, G'zz = 1, C s xy — —Cy X — is/2, and 
otherwise. 

Each of the three terms in Eq. (|A1I) is a convolution. 
By applying the convolution theorem, we obtain 

Ip(r,<p,K) = y2f- 1 \K a p(k)fM a (T)} , (A5) 

l J z—0 

a 

where T stands for the Fourier transform. The functions 
K a p{k) — TDapir) can be written in cylindrical coordi- 
nates (k, 8, k z ) as 

2 

K a0 (k,6,k z ) = nY, K e a p(kz/k)e m , (A6) 
e=-2 

where K% = -K± 2 = ±iK^ = L(0 = 1/(1 + £ 2 ), 
Kg = ±iK^ = 2£L(0, K° xx - K% y = -4/3 + 2L(£), 
K° zz = -4/3 + 4£ 2 £(0, K l ap = K* a , and otherwise. _ 
The Fourier transforms of the components of magneti- 
zation in Eq. (| A5|) are readily evaluated with the aid of 
the Jacobi- Anger expansion, and the result is 

FM a {v) = (2^) 2 ]T C s aj (-i) n H^(k)e ine 5( K (s+n)-k z ), 

7.s,n 

(A7) 

where H™(k) = J rdrJ n (kr)g™(r) is the nth order Han- 
kel transformation of gUj. By substituting Eqs. (|A6[) and 
(|A7p into Eq. (|A5|) and carrying out the integration over 
k z , we obtain 

I (r, ( p,K) = 2n 2 J2 CS a -y(- i ) n ( A8 ) 



7 



where q = n(s+n) / k and J-^d denotes the inverse Fourier 
transform in the z — plane. Again, the inverse Fourier 
transforms can be evaluated by applying the Jacobi- 
Anger expansion, resulting in 



/ z "' ±1 =2(ff x "T^") q±L(q ± ), 



fl 



2h: 



2 9o 2 L((Zo) 



(A13) 
(A14) 



(A9) 



where h^(r) = J kdkJ n+ i(kr)f^ i (k, n) stands for the 
inverse Hankel transformation of order n + 1 of the func- 
tion ff(k,K) = L !7 /;^(#^)- By denoting 
q± = k{u zb l)/k and q$ = Kn/k, these functions can be 
written as 



fn,±2 
J x 

J X 



TiH^)L{q ± )=±if, ! 



= 2H?q L(q a ) = ±if y 



Tt,±2 

y ' 



(A10) 
(All) 



±i/"> u ,(A12) 



where /™ : ° 
Eq. (|An|) . 



and 



fn.O 



are given implicitly through 



In the case of an axially polarized cylindrically sym- 
metric state, ipk = \/n(r)5ik, the GP equation, Eq. ([1}, is 
reduced to a scalar equation. In this simple case, g® yields 
the only non- vanishing term in Eq. (IA3|) . which, together 
with the fact that q = 0, implies that /°' ° = -4ff°/ 3 
is the only non- vanishing function in Eqs. (IA10|) — (| A14|) . 
Substitution into Eq. (|A9I) shows that Eq. (QJ is reduced 
into the scalar equation if we replace g — > g + g s — -j-<?dj 
which is readily proved also by direct evaluation of 
Eq. {Mi- 
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